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ABSTRACT 

We present new Arecibo 21-cm line measurements of the temperatures of interstellar gas. Our tem- 
peratures for the Cold Neutral Medium (CNM) are significantly lower than previous single-dish results 
and in very good accord with theoretical models. For warm gas at T > 500 K, we find a significant 
fraction of gas (> 47%) to lie in the thermally unstable region 500 — > 5000 K; moreover, about 60% of 
all the neutral atomic gas has T > 500 K. Large amounts of thermally unstable gas are not allowed in 
theoretical models of the global interstellar medium. 

Subject headings: ISM:atoms — ISMrgeneral — radio lines: ISM 
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1. INTRODUCTION 

In 1977 McKee & Ostriker (MO) extended the two-phase 
model of Field, Goldsmith, & Habing (1969) by including 
supernovae (SN) . Their interstellar medium (ISM) has four 
phases, if one includes the Warm Ionized Medium, and is 
dominated by individual SN explosions. The Hot Ionized 
Medium (HIM) fills the interior of SN remnants and pow- 
ers their blast waves, which sweep up the gas inside the 
bubbles and pile it into the shells. Soon this shocked gas 
starts to cool and recombine rapidly, forming the Cold 
Neutral Medium (CNM). Soft X-rays produced by adja- 
cent HIM penetrate the outsides of CNM clouds, heating 
the gas to form the Warm Neutral Medium (WNM). 

The swept-up gas consists of two distinct neutral phases 
in physical contact, each in thermally stable equilibrium 
with equal thermal pressures. It is difficult to study such 
gas in individual SN remnants, but for superbubbles the 
general picture of hot gas inside a swept-out volume sur- 
rounded by dense walls is well corroborated by modern 
multiwavelength studies, for example in Eridanus (Heiles, 
Haffner, & Reynolds 1999). 

Theoretical CNM and WNM temperatures are derived 
by calculating the equilibrium temperature as a function 
of thermal pressure. There exist two stable ranges of equi- 
librium, the CNM and the WNM with temperatures ~ 50 
and ~ 8000 K, separated by a region of unstable tempera- 
tures (Wolfirc et al 1995; WHMTB). Stable thermal pres- 
sure equilibrium, and thermal pressure equality between 
the phases, is a cornerstone of the MO theory. 

Observational temperatures are derived in different ways 
for the CNM and WNM. For the CNM, the most accurate 
temperatures come from comparing 21-cm line absorption 
and emission profiles. The line opacity oc ^; only the 
CNM produces significant absorption lines, while the CNM 
and WNM both contribute to emission. In contrast, most 
WNM temperatures are derived from linewidths. 21-cm 
linewidths provide upper limits to temperature, while in 
some cases combining them with heavy element linewidths 
from optical/UV spectra provides actual temperatures (see 
review by Heiles 2000). 

Kulkarni & Heiles (1987; KH) summarized tempera- 
ture measurements of the CNM and WNM, confining their 



analysis to 21-cm line absorption/emission measurements 
because very few UV data were available at that time. 
For the CNM, they discussed single-dish results in terms 
of the conventional interpretation of clouds within which 
the temperature increases outwards. As briefly discussed 
below, each derived temperature is the lowest one in its 
cloud. Expressed as histograms, these coldest-cloud de- 
rived temperatures are broadly distributed over the range 
20 ^ 300 K (Mebold et al 1982, MWKG; Dickey, Salpeter, 
& Terzian 1978, DST; Payne, Salpeter, & Terzian 1982, 
PST). In addition, there is a weakly significant statistical 
relationship between the derived CNM temperature and 
the 21-cm line opacity, called the T — t relation. For the 
WNM, KH found that the limited data supported a lower 
limit ^ 5000 K, but they cautioned that the result needed 
confirmation. 

In contrast, interferometric maps of a field around 
3C147 (Kalberla, Schwarz, & Goss 1985, KSG) provide 
a completely different picture. For the CNM, the compo- 
nents are roughly isothermal and have colder temperatures 
(34 ^ 74 K). For the WNM, which contributes - 80% of 
the mass, temperatures lie in the thermally unstable range 
(500 — > 2000 K). The high angular resolution should make 
these results reliable. However, such results are available 
in only a few fields. The fact that they conflict with the 
single-dish results is disturbing. The single-dish results are 
the ones always quoted in reviews because of their much 
larger statistical sample for the CNM; it is worth noting, 
though, that for the WNM the single-dish sample is no 
larger than the interferometric one. 

Here we present a statistical summary of new single- 
dish temperatures. For the CNM they are derived from 
HI absorption/emission line data; for the WNM they are 
upper limits based on linewidths. We introduce a new 
analysis technique for absorption/emission observations. 
Our results are consistent with the 3C147 results. For the 
CNM, we find lower temperatures than previous single- 
dish workers; these are more in line with theoretical pre- 
diction for the CNM. For the WNM we find that most of 
the gas WNM is at thermally unstable temperatures; this 
disagrees with the MO theory. 

2. NEW 21-CM LINE OBSERVATIONS 
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2.1. Observations and reduction technique 

Heiles and Troland (2001) are performing a survey of 
Zeeman splitting of HI absorption lines with the Arecibo 
telescope^. These data have long integration times, which 
produces excellent signal/noise and makes them unsur- 
passed for obtaining temperatures. Our results arc more 
accurate and cover more sources than the best previous 
surveys, which are the Arecibo work by DST^ and PST, 
and the Bonn/NRAO work by MWKG. Here we report 
on 24 sightlines, 19 of which have |6| > 20° but otherwise 
arc randomly selected within Arecibo's declination range 
~ 0° ^ 39°. 

Each absorption spectrum consists of very obvious ve- 
locity components and we represent their optical depths 
by a set of N Gaussians. Thus we least squares fit the 
observed spectrum = e~^, where 



^0 



-1_ p-[{V-Von)/SVnr 



(1) 

here is the absorption profile divided by the contin- 
uum source strength and (Vbn, 5Vn) are (central velocity, ^ 
width) of component n. We assume that each component 
is an independent physical entity, and is isothermal This 
is consistent with the findings of KSG. 

Each HI emission spectrum contains structure but, also, 
is wider than its associated absorption spectrum, as is 
well known. We assume that the emission spectrum 
Tern = TcNM + TwNM, where Tcnm is the contribution 
from the aforementioned CNM components. Twnm is the 
contribution from K additional wide Gaussians to repre- 
sent the WNM; K is small number and often just one. The 
spin temperatures in these additional components are so 
high that they have negligible optical depth and produce 
no easily discernible features in the absorption spectrum. 

In least-squares fitting the emission, wc include the ab- 
sorption of more distant CNM Gaussians by less distant 
ones. Letting T„ be the spin temperature of component 
n, which is also the kinetic temperature, 



CNM 



= S^-iT„(l - e- 



(2) 



where the subscript m with its associated optical depth 
profile Tm represents each of the M CNM clouds that lie 
in front of cloud n. For multiple absorption components, 
we experiment with all possible orders along the line of 
sight and choose the one that yields the smallest residuals. 
We also include the absorption of each WNM component 
by the CNM by assuming that a fraction J^k lies in front 
of all the CNM and is unabsorbed, with the rest all lying 
behind; thus 

Twnm = S^-i[J-fe + (l- J-fc)e-^]rofee-[(^-^°'=/*^'=l' , (3) 
where the subscript k represents each WNM component. 
Note that Tok is a brightness temperature, not a kinetic 
temperature. In most cases J^k is indeterminant and we 
can only distinguish between the two extremes J-'k = (0, 1). 
The differences for different orderings are sometimes not 
statistically significant but nevertheless lead to differences 
in the derived CNM temperatures. These differences re- 
flect the uncertainties in r„ more than the conventional 

1 The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell University under a 
cooperative agreement with the National Science Foundation. 

^ We find about half the DST absorption profiles to exhibit large differences from ours, with DST components being wider and multiply-peaJjed; 
some of DST's profiles were corrupted by local oscillator stability problems (Dickey, private communication). This completely explains the 
disagreement of Greisen and Liszt (1986) with DST for 3C348. 



errors derived from least squares fits. Additional uncer- 
tainties can occur if there are unresolved subcomponents. 
We defer discussion of these details to the more compre- 
hensive paper (Heiles and Troland 2001). 

Previous single-dish authors, in contrast, implicitly as- 
sume that clouds are not isothermal. They derive the spin 
temperature of a cloud at the peak of its absorption pro- 
file. Thus each point on their histograms represents the 
lowest derived temperature for that particular cloud. This 
temperature, however, is not the coldest temperature in 
the cloud, because the line of sight also passes through 
warmer gas. 

2.2. Sample result: 3C18 

Figure 1 exhibits the results for 3C18 [located at (^, h) = 
(119°, —53°)], which is a simple profile and good for an 
illustrative example. In the top panel the solid line is 
the observed absorption spectrum which we fit with 
the three CNM components whose depths and halfwidths 
arc indicated; the dash-dot line is the fit. In the bottom 
panel the solid lino is the observed emission spectrum T^m ■ 
The dashed curve is Tcnm', the dotted is T\\-^\m fit with 
K = 1, which is unabsorbed by the CNM because the low- 
est residuals are obtained with = 1. The full fitted curve 
is the sum, shown as dash-dot, which is a good fit except 
in the extreme line wings where stray radiation makes the 
data suspect (e.g. Hartmann & Burton 1997). 

For 3C18, the WNM component has halfwidth 10.0 km 
s~^, which corresponds to purely thermal broadening at 
T = 2200 K; this is an upper limit on the kinetic tempera- 
ture Tk- For the three CNM components, Icft-to- right on 
Figure 1, the halfwidths are (2.5±0.03, 6.3±0.07, 1.1±0.15) 
km s~^ and spin temperatures are 32±1, 43±6, and 46±9 
K. The ratios of total linewidth to thermal linewidth are 
(2.04 ± 0.07, 4.50 ± 0.63, 0.75 ± 0.18); the Icr uncertainty 
on the last ratio is statistically consistent with a ratio 
> 1, as it must be. The (WNM, CNM) components con- 
tribute NiHI) = (3.2, 1.8) X 10^0 cm^2^ respectively. The 
WNM/CNM ratio is ~ 1.8, which is close to our global 
average. 



2.3. Our ensemble of WNM temperatures 

Our 49 WNM components have linewidths that corre- 
spond to upper limits on the kinetic temperature. The top 
two panels of Figure 2 exhibit histograms of these limits, 
one for number of components and one for column den- 
sity. Not included on these histograms is one absorption 
component for which the spin temperature was derived: 
T,p,„ = 725 K, N{HI) = 1.4 x 10^° cm^^. Including this, 
20 of the WNM components (40%) have Tk = 500 ^ 5000 
K. These contain > 47% of the total WNM column den- 
sity; this is a lower limit because the WNM temperatures 
are upper limits. Because these components are not visi- 
ble in absorption, their spin temperatures exceed ~ 500 K. 
This range, 500 5000 K, is approximately the thermally 
unstable range that separates CNM from WNM. 
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2.4. Our ensemble of CNM temperatures 

Our 86 CNM temperatures are derived from absorp- 
tion/emission data and arc values, not upper or lower lim- 
its. The bottom two panels of Figure 2 exhibit histograms 
of these temperatures. The two histograms exhibit broad 
peaks in the range T = 25 ^ 75 K. 47 of the CNM compo- 
nents (54%) have temperatures in this range; these contain 
61% of the total CNM column density. We also see colder 
gas: 10 components (11%) containing 5% of the mass have 
T = 10 25 K. We discount the four sm3l\-N{HI) com- 
ponents having T < 10 K: they are weak and the tempera- 
tures have large errors. We find no support for the weakly 
significant T — t relation reviewed by KH. 

2.5. Ratio of CNM and WNM components and mass 

For WNM gas (T > 500 K), we found 49 components 
with total N{HI) = 107 x 10^° cm'^. For CNM gas 
(T < 200 K) we found 80 components having a total 
N{HI) = 65 X 10^0 cm~2 ^j^^, mildly ambi guous range 
T = 200 — + 500 K we found 6 components with total 
N(HI) = 7.5 X 10^0 cm-2. Thus, the ratio of CNM to 
WNM is, in terms of number of components, 1.6; in terms 
of mass, 0.60. Overall, our results indicate that about 60% 
of all the neutral atomic ISM is WNM, with T > 500 K. 

2.6. Temperatures from optical/UV absorption line 
observations 

To derive kinetic temperatures from an atomic opti- 
cal/UV absorption line, one decomposes it into Gaussian 

components. Then one docs the same with the 21-cm emis- 
sion line toward the star, fixing the central velocities to be 
the same. The comparison of line widths separates the 
thermal and turbulent broadening. The derived temper- 
atures are upper limits because the HI line comes from a 
much larger angular area so the nonthermal component of 
its width may be larger than that of the heavy element 
lines. Such temperatures are probably the best one can do 
for the WNM, but for CNM gas they are much less accu- 
rate than those derived from 21-cm absorption/emission 
line data. 

Spitzer & Fitzpatrick (1995) and Fitzpatrick & Spitzer 
(1997) use this technique towards two high-latitude stars 
and derive temperatures and column densities for 21 dif- 
fuse neutral components. Of this total, 3 components have 
T > 5000 K, 13 have T < 500 K, and 5 have T in the un- 
stable 500 — > 5000 K range; thus, 24% of the components 
are thermally unstable; these contain 63% of the mass. 

One can derive excitation temperatures of the low-J 
states of H2 using UV absorption lines (ShuU et al 2000; 
Spitzer, Cochran, & Hirshfeld 1974). These tend to agree 
with previous 21-cm line temperatures and are systemat- 
ically higher than ours. If our temperatures are correct, 
then the low-J states are nonthermally populated, as are 
the high-J ones. 

3. DISCUSSION AND COMPARISON WITH THEORY 

3.1. The WNM 

Both our new HI and the optical/UV observations show 
that much of the WNM — at least 47% — lies at temper- 
atures that are unstable to isobaric perturbations. Our 
Arecibo data show this departure from thermal stability 



in a statistically convincing manner. Previous 21-cm line 
studies have hinted at this result. In emission/absorption 
studies, MWKG decomposed emission line profiles into 
Gaussians, with similar results; however, they didn't ex- 
plicitly point out this departure. Vcrschuur & Magnani 
(1994) and Heiles (1989) analyzed emission profiles and 
found numerous components with widths in this range, 
but without absorption data could not conclusively state 
that the kinetic temperatures were indeed so high. 

The large fraction of WNM in the thermally unstable 
regime violates a fundamental cornerstone of equilibrium 
ISM models such as MO, which all rely on thermal pressure 
equilibrium to push the gas into one of the thermally stable 
CNM or WNM phases. This result seems to push us to- 
wards other types of model. Two possibilities include time- 
dependent models such as the supernova-dominant model 
of Gerola, Kafatos, and McCray (1974) and turbulence- 
dominated models such as Vazquez-Semadeni, Gazol, & 
Scalo (2000). 

3.2. The CNM 

Figure 2 exhibits the histogram of derived spin temper- 
atures for all CNM components. Both most of the compo- 
nents and most of the mass have T = 25 ^ 75 K. This is 
in marked contrast to previous results, where histograms 
were broad over the ranges 20 — > 140 K (MWKG) and 
50 ^ 300 K (DST, PST). Our range is narrower and, 
moreover, temperatures extend to very low values, with 
significant contributions down to T = 10 K. 

The peak above T ~ 25 K agrees very well with the 
high angular resolution results of KSG and, also, theory. 
WHMTB included all known processes in calculating their 
standard model, for which the CNM equilibrium tempera- 
tures range from 25 — > 200 K (the corresponding densities 
are uhi ^ 1000 — » 4 cm~^). Our observed temperature 
range is smaller and corresponds to nHi ^ 250 — > 20 cm~^ 
and f = 10000 1500 cm ^ K. These numbers are in ac- 
cord with ISM pressure measurements (Jenkins, Jura, & 
Lowenstein 1983). 

Temperatures as low as our 10 — > 20 K range can oc- 
cur in the absence of the PAH-type grains that produce 
grain heating (WHMTB; Bakes and Tielens 1994). In this 
case, heating is by photoionization of Carbon and cooling 
by electron recombination onto ionized Carbon (Spitzer 
1978). Such cold (and even colder) gas was invoked by 
Heiles (1997) to help understand the existence of tiny-scale 
atomic structure; the present results are encouraging for 
that interpretation. 

It is a pleasure to acknowledge conversations with John 
Dickey, Alex Lazarian, Jeff Linsky, Chris McKee, Paolo 
Padoan, Tom Troland, Enrique Vazquez-Semadeni, and 
Ellen Zweibel. This work was partly supported by NSF 
grant AST9530590 to the author. 
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Fig. 1. — 21-cm line absorption (top) and emission (bottom) spectra for 3C18. In the absorption spectrum, the solid line is data and the 
dash-dot line the fit; crosses indicate the central e~'^'s and halfwidths of the three Gaussian component parameters. In the emission spectrum, 
the solid line is the data; the dashed line is the contribution from the three CNM components, the dotted line from the WNM component, 
and the dash-dot line their sum. Dash-dot lines in both figures are the fits and are so close to the data that they axe hard to distinguish. 
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Fig. 2. — Histograms of derived temperatures. The top two frames axe upper limits on kinetic temperature for the WNM derived from 
line widths; of these, the top gives the number of Gaussian components and the bottom the HI coluiriii density in units of 10^'^ cm~^. Off 
the graphs to the right, with upper limits exceeding 10* K, are 21 components (44%) containing N{Hr) = 45 X 10^° cm-^ (43%). The 
bottom two frames are values (not limits) of spin temperature derived from absorption/emission data. Off the graphs to the right, with spin 
temperatures exceeding 200 K, axe 8 components (9%)containing N{HI) = 37 X lO^" cm-^ (36%). 



